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Abstract 

We show that in tilting perfect fluid cosmological models with an 
ultra-radiative equation of state, generically the tilt becomes extreme 
at late times and, as the tilt instability sets in, observers moving with 
the tilting fluid will experience singular behaviour in which infinite 
expansion is reached within a finite proper time, similar to that of 
phantom cosmology (but without the need for exotic forms of matter). 

PACS: 98.80.Jk, 04.40.Nr 

Cosmological data, including galaxy, CMB and supernovae observations 
PP, seem to be consistent with a cosmological constant or a dark phantom 
energy with effective equation of state parameter 7 < 0 [2]. We shall study 
cosmological models with a tilting, but otherwise conventional, perfect fluid, 
and show that the models have dynamical behaviour similar to that of phan¬ 
tom cosmology, but without the need for any exotic forms of matter and 
consequently avoiding the pathologies in these models, such as the existence 
of ghosts. 

For spatially homogeneous (SH) Bianchi cosmologies models, the universe 
is foliated into space-like hypersurfaces [3], and there are two naturally de- 
hned time-like vectors: the unit vector held, ? 7 ,“, normal to the group orbits, 
and the four-velocity, of the perfect huid. If is not aligned with n“, the 
model is called tilted (and non-tilted or orthogonal otherwise) |1]. Usually, 
the kinematical quantities associated with the normal congruence ?7,“ of the 
spatial symmetry surfaces, rather than the huid how are used as vari¬ 
ables. Following jl] , a tilt variable v is introduced, so that in an orthonormal 
frame where = (1, 0, 0, 0), we have 


\/l — v'^ 


(l,u,0,0) 


( 1 ) 
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We will assume a perfect-fluid matter source with p = (7 — 1 )/! as equation 
of state, where fi is the energy density, p is the pressure, and 7 is a constant. 
Causality then requires 7 to be in the interval 0 < 7 < 2 . A positive 
cosmological constant may also be included in the models. Such SH tilted 
cosmologies with a 7 -law perfect fluid source have been studied by a number 
of authors ini[7iEiEi[iniiiiiE|. It is known that the tilt can become extreme 

^ 1) asymptotically to the future. In particular, let T be the proper 
time as measured along the fluid congruence, and let us dehne the quantity 
AT: 

AT = / —"s/l — v^dr (2) 

J TO H 

where r is a dynamical time variable dehned in terms of the clock time t and 
the Hubble scalar H hy ^ = H. If AT is hnite, then the fluid congruence 
is future incomplete: the fluid observers will reach inhnite expansion within 
hnite proper time. Therefore, in spite of the fact that such a spacetime 
may be future geodesically complete ^ 21 ; the worldlines dehned by the huid 
congruence u°‘ may become null with respect to the normal congruence n“, 
sometimes so quickly that this occurs within hnite huid proper time. We note 
that for a SH cosmology for which 1 (which is a necessary condition 

for AT to be hnite) the Hubble parameter of the huid congruence, if fluid also 
diverges. 

The locally rotational symmetric (LRS) Bianchi type V perfect huid mod¬ 
els were studied in j7j. The global structure of Bianchi type V models was 
studied in more detail in jB]. It was shown that for models with 4/3 < 7 < 2, 
to the future the tilt can become extreme (in a hnite time as measured along 
the huid congruence) and, as T approaches a hnite limiting value, the huid 
worldlines become null with respect to and its length scale i and some of 
its kinematic variables (i.e., the expansion, the shear, the acceleration, and 
the vorticity) diverge, but its matter density and curvature scalars tend to 
zero. Although this peculiarity only occurs for ultra-radiative huids, which 
are unlikely to dominate the late universe, we will show that this behaviour 
can even occur for universes where ultra-radiative huid is not dominant. 

We shall hrst show that the huid proper time is hnite as the solution 
approaches its asymptotic state (the extremely tilted Milne equilibrium point 
M“, dehned by E = 0 , A = 1 , n = — 1 , H = 0 , Ha = 0 , and with deceleration 
parameter O' = 0) for 4 / 3 < 7 < 2 in the absence of a cosmological constant 
[ 12 ]. The decay rates for the Hubble scalar H = 9/3 and the quantity a/I — 
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are H (x e , \/l — oc exp [—(57 — 6)r /(2 — 7)], whence 


AT cx 



■00 


exp 



( 3 ) 


which is finite for 4/3 < 7 < 2. 

Spatially homogeneous cosmological models with a positive cosmological 
constant were investigated using dynamical systems methods in na. extend¬ 
ing the tilted LRS Bianchi type V analysis of [Zj to the A 7 ^ 0 case. A de 
Sitter point with extreme tilt is the future attractor for 7 > 4/3. Therefore, 
in general, for Bianchi type V models with 7 > 4/3, the tilt again becomes 
extreme at late times and the fluid motion is no longer orthogonal to the 
surfaces of homogeneity. 

Although it is known that expanding non-type-IX Bianchi models with a 
positive cosmological constant isotropize to the future (cosmic no hair theo¬ 
rem) [ 13 ], and that this result applies to tilted models, the isotropization is 
with respect to the congruence normal to the homogeneous symmetry sur¬ 
faces - not the fluid congruence |T3]. Thus in the Bianchi type V models the 
space-time generically becomes de Sitter-like, in accordance with the cosmic 
no hair theorem, but since the tilt does not die away, isotropization of the 
cosmology does not occur with respect to the fluid congruence. 

This is a generic feature of spatially homogeneous models. In general, SH 
models are not asymptotically isotropic. The spatially homogeneous models 
do isotropize to the future in the presence of a positive cosmological con¬ 
stant m By investigating the asymptotic behaviour of a SH model with 
a cosmological constant and a tilted perfect fluid with p = (7 — l)/i as the 
future de Sitter model (with an extremely tilted perfect fluid) attractor is 
approached we obtain H x Hq and: AT oc exp [—(37 — 4)r/(2 — 7 )] dr , 
which is finite for 4 / 3 < 7 < 2 . In the absence of a cosmological con¬ 
stant, in general SH models are not asymptotically isotropic. Neverthe¬ 
less, these models can spend a long time close to a fiat Friedmann model 
corresponding to a saddle point. Considering one non-tilting perfect fluid 
with p± = (j± — l)ju± and one tilting fluid with p = (7 — l)/i (with 7 ^ < 
7 ), as the Friedmann equilibrium point is approached we obtain: AT oc 

exp [—(37 — 4 — |7±(2 — 7 ))r/(2 — 7 )] dr. This means that for solutions 
spending a finite (but arbitrarily long) time close to the saddle T, exhibiting 
quasi-isotropic behaviour consistent with observations, the Hubble parameter 
as measured by the fluid, ffauid, can become arbitrarily large. 
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There will be different anisotropic asymptotic end-states depending on 
Bianchi types. In the Bianchi type VIIq model with a tilted 7 -law perfect fluid 
mcni, the future asymptotic state for 7 > 4/3 was found to be anisotropic 
and extremely tilted with H oc and AT oc exp [(8 — 57 )r/(2 — 7 )] dr. 
For models with 7 > 8/5, we find that this integral is finite. The Bianchi type 
VIII models are asymptotically extremely tilted for 1 < 7 < 2 and the asymp¬ 
totic solution is an extremely Weyl-curvature dominated model Ha CHI, with 
H oc r5exp(—|r), and AT oc r 4 ( 2 - 7 ) exp [— 3(87 — 4)r/2(2 — 7 )] dr, 
which is hnite for 4/3 < 7 < 2. 

It is known that for general spatially inhomogeneous perfect fluid models 
with a cosmological constant, the de Sitter solution with extreme tilt (where 
the tilt refers to the fluid tilt with respect to a congruence with an acceleration 
that tends to zero) is locally stable for 4/3 < 7 < 2 [1211121 • From Equations 
(3.43) and (3.28) of [ 12 ], we have that H oc Hq and 

AToc^ 

which is finite for 4/3 < 7 < 2, as required. We note that for 4/3 < 7 < 2, the 
generic behaviour is n —>■ 1 , which implies that inflation does not isotropize 
an ultra-radiative fluid. Even if inflation is turned off after a certain number 
of e-foldings, FTfluid can become arbitrary large. Nevertheless, as noted above, 
this result does not contradict the cosmic no-hair theorem dl 

Therefore, we have found that, for 7 > 4/3, the fluid congruence becomes 
null with respect to the normal congruence in finite fluid proper time, and a 
’kinematic singularity’ develops for the fluid congruence. To fully understand 
the behaviour of these models and their physical properties, the dynamics 
need to be studied using a formulation adapted to the fluid (i.e., utilizing 
a fluid-comoving frame). By using the boost formulae relating the normal 
and fluid congurences this singular behaviour for the fluid congruence can be 
confirmed. 

This mathematical instability might lead to some interesting physics. Ex¬ 
panding universes that come to a violent end after a finite proper time have 
arisen in a different context EHHOj. Models with a constant equation of state 
parameter 7 < 0, dubbed “phantom energy”, lead to a singularity commonly 
called the big rip. In this paradigm, during the cosmic evolution the scale 
factor grows more rapidly than the Hubble distance and consequently blows 
up in a finite proper time, and is typically characterized by a divergent pres- 
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sure and acceleration. As the big rip singularity is approached, both the 
strong and weak energy conditions are violated. For the so-called sudden 
future singularities [201 the strong energy condition needs not be violated, 
yet a future singularity forms within hnite time. In our examples, no energy 
conditions are violated; in fact, the energy density in the fluid frame, /ifiuid, 
tends to zero as the singularity is approached. 

The Hubble scalar for the fluid frame can be computed from the boost 
formula for H, and is of the form iiffluid = BH/\/l — v‘^ j2I]. The above 
examples lead to a diverging iffluid, with i^guid ot = 

Hr in the limit. We then examine the value of the deceleration parameter 
gfiuid and compare with the critical value of —1. Equivalently, we can compare 
with a phantom fluid in a flat, isotropic model using the effective equation 
of state parameter 7 eff, given by 


beff 3 *?fluid)) 


(5) 


with corresponding critical value of 0 . 

In the general case of inhomogeneous cosmological models with a cosmo¬ 
logical constant, the de Sitter model with extreme tilt is the future asymptotic 
state for 7 > 4/3, and we have 


3 

muid = “ 2(37 - 4) - 1 < -1. 


( 6 ) 


Consequently, as the de Sitter asymptotic state is approached the dynami¬ 
cal effect of the ultra-radiative perfect fluid congurence behaves similarly to 
a phantom energy in an isotropic and spatially flat spacetime. Note that 
inflation does not stop the big rip from occurring. 

Let us now briefly consider the SH models discussed above. For the 
LRS Bianchi type V models, on the approach to the extremely tilted Milne 
solution (for 7 >4/3) we hnd that 


3 

gfluid = - 4) - 1 < - 1 . 


For the Bianchi type VIII model with 7 > 4/3 we again obtain 
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For the Bianchi type VIIq model with 7 > 8/5 we obtain 

3 

<?fluid = --(57 - 8 ) - 1 < - 1 ; 

this difference arises dne to the fact that the type VIIq models are geometri¬ 
cally more special than Bianchi type VIII models. Similarly, for the two-ffnid 
example, we have 

*?fluid ^ 1) 

for 7 > ^ similar behavionr also occnrs for Bianchi type VII/j models 

(see |2I]). In the generic SH models (snch as, for example, the Bianchi type 
VIII model) the nltra-radiative perfect ffnid effectively behaves dynamically 
like a phantom energy in the sense that the length scale and the Hnbble scalar 
diverges as the fntnre asymptotic state is approached. In the more special 
examples, the reqnirement is that the threshold eqnation of state mnst be 
eqnal to or higher than that of radiation. 

Therefore, as the fntnre asymptotic state is approached the nltra-radiative 
perfect ffnid effectively behaves like a phantom energy in an isotropic and 
spatially flat spacetime. It is important to note that the energy conditions 
of the perfect ffnid are nowhere violated. 

Let ns discnss the physical conseqnences of this dynamical behavionr in 
a little more detail. One can consider models that spend a period close to 
isotropy (i.e., close to a Friedmann saddle point), with a small tilt. There¬ 
after, the models begin to evolve away from isotropy. Since the tilt is non¬ 
zero, for 7 > 4/3 the models generically evolve towards an asymptotic state 
with extreme tilt. As the tilt instability sets in dnring the transient regime, 
observers moving with the tilting ffnid will experience a transition from a 
decelerating expansion to an accelerating expansion, and later, extremely 
accelerating expansion mimicking that of a phantom cosmology. Moreover, 
nnlike in conventional phantom cosmology, in the models stndied here there 
is no need for any exotic forms of matter; conventional matter which is tilt¬ 
ing snfhces. In a braneworld approach, accelerating nniverses can also resnlt 
withont a cosmological constant or other form of dark energy E3- Indeed, 
other pathologies, snch as the existence of ghosts, are avoided in the mod¬ 
els described here. This is also the case in alternative models to phantom 
cosmology which resnlt from alternative theories of gravity, theories with non- 
minimal conplings, and models in which the dark energy and qnintessence 
field interact [221 • In addition, as noted above, dne to the existence of fn- 
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ture attractors with extreme tilt the dynamical behaviour described here is 
generic. 

Although the dynamics in the comoving dark energy models and the tilt¬ 
ing models are qualitatively similar, the actual quantitative (physical) pre¬ 
dictions of the two different models may differ (due to the different physical 
transient time scales in each model). For example, the age of the universe 
as measured in the two models may differ. As an example, let us calculate 
numerically the age in a LRS Bianchi type V model which starts near a flat 
FL model at the time of decoupling. In a conventional comoving dust model 
(with no dark energy) with present energy density fIdust,o ~ 0 .2, the maxi¬ 
mum age is approximately 11.8 billion years (which, as is well known, is on the 
low side compared with the ages of the oldest astrophysical objects in the Uni¬ 
verse). In a comoving dust model with present energy density fldust,o ~ 0.2 
and a cosmological constant Ua,o ~ 0.8, the estimated age is about 15.0 bil¬ 
lion years, which is about a 27% increase over the age for the dust model and 
consistent with current observational data. In a LRS Bianchi type V model 
with a comoving dust model with fldust,o ~ 0.2 and a (second) tilting perfect 
fluid with a total effective energy density fltiit.o ~ 0.8, then according to the 
tilted fluid observer with present energy density Udust.o,fluid frame = 0.2, the 
maximum age as measured in the tilting fluid frame appears to be (depending 
on the initial conditions in the various numerical experiments) about 13.2 bil¬ 
lion years (this particular value occurs for the initial conditions: Aq = 0.641, 
S+_o = 0.0056, V = —0.443; Uq = 0.0111). This is a 12% increase over the 
original age of 11.8 billion years (and marginally consistent with observa¬ 
tions). The situation is expected to be similar in the more general Bianchi 
VIII models, although the quantitive predictions will depend on the precise 
initial conditions and it is possible more hne tuning may be necessary. 

Therefore, although qualitatively the calculated age of the universe in 
both the dark energy models and the tilting fluid models are similar in that 
the age is greater than in the conventional models (without dark energy or a 
second tilting fluid), quantitatively the increased age appears to be greater in 
the dark energy models. However, the tilted fluid models are still physically 
viable. In future work we shall further study physical predictions of the 
tilting fluid models. In addition to the age problem, it is also of interest 
to investigate the effect of tilt on cosmic microwave background radiation 
observations and whether these models offer a possible explanation for the 
various anomalies on large angular scales found in the WMAP data EH] . 
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